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Abstract 

The Dirac equation in (1+1) dimensions with a non-local VT- 
symmetric potential of separable type is studied by means of the Green 
function method: properties of bound and scattering states are derived 
in full detail and numerical results are shown for a potential kernel of 
Yamaguchi type, inspired by the treatment of low-energy nucleon- 
nucleon interaction. 
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1 Introduction 

Since the pioneering papers by Bender and coworkers [T], [2], the study of 
non- Hermit ian Hamiltonians invariant under space-time reflection has devel- 
oped into a branch of quantum mechanics in its own, called "PT-symmetric 
quantum mechanics. The large majority of analyses have been devoted to 
bound state problems, where the observation that "PT-symmetric Hamiltoni- 
ans with eigenfunctions that are eigenstates of VT have real spectra has led 
to Hermitian-equivalent formulations, where one can define a linear operator 



C, commuting with Hamiltonian H and with VT , that permits construct- 
ing time-independent inner products with positive-definite norms of the form 
/_ °° ^ i x ) CVTty (x) dx (see Ref. [3] for a review). 

'PT-symmetric quantum mechanics has a close connection with the more 
general quasi-Hermitian quantum mechanics [I], [5], [6], where H is called quasi- 
Hermitian if it satisfies the intertwining relation = r] + Hr]^ 1 , with r/ + a 
positive- definite Hermitian operator called the metric operator, playing a role 
analogous to CV. 

While bound states of PT-symmetric Hamiltonians are nowadays well 
understood, many more questions remain open in the treatment of scatter- 
ing states: for instance, it has been shown in Ref. [7] that even simple local 
potentials require introducing non-local metric operators and non-standard 
boundary conditions with progressive and regressive waves not only in the 
initial but also in the final state. Even if the latter feature might be re- 
moved by an appropriate choice of the metric operator, called quasi-local 
in Refs. j8],[9], a satisfactory general approach has not been formulated yet. 
This is why the majority of studies on scattering by PT-symmetric Hamilto- 
nians has been made within the framework of standard quantum mechanics, 
breaking unitarity of the corresponding scattering matrices (see Ref. [11] and 
references therein). Even at this effective level, PT-symmetric potentials are 
peculiar, in the sense that, depending on their parameters, they may behave 
as absorptive for progressive waves and generative for regressive ones (or 
viceversa), a property called handedness in Ref. [10], or it may happen that 
they are neither absorptive, nor generative, because the sum of the square 
moduli of transmission and reflection coefficients may be smaller than one, 
or greater than one in different intervals of incident energy; they can even 
conserve unitarity when the asymptotic wave functions are eigenstates of 
VT: in this latter case they are necessarily reflectionless[TTJ. As is known, 
the reflection of progressive (left-to-right) and regressive (right-to- left) waves 
is quite asymmetric (Rl-,r 7^ Rr^l) already in the case of local potentials, 
where the transmission is the same; in the case of non-local potentials[12j, the 
transmission is asymmetric, too (T L _> R 7^ T R ^ L ). Indeed, non-local poten- 
tials have more subtle PT-transformation properties than local potentials, 
for which T- invariance and Hermiticity requirements coincidepTj. 

The scenario is even richer in relativistic quantum mechanics, where, 
again, the majority of studies have been dedicated to bound states of VT- 
symmetric potentials in the Klein-Gordon and Dirac equations in (1+1) 
space-time dimensions. Limiting ourselves to the Dirac equation, of interest 
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to the present work, we may quote the pseudo-supersymmetric description [13], [2] 
of scalar or pseudo-scalar local potentials with exact, or spontaneously bro- 
ken PT symmetry, the PT-symmetric version of the generalized Hulthen 
vector potential[15j , the combinations of scalar (position-dependent mass) 
and vector potentials of Refs. [IS] . |17j . [IB] . [19] . 

Making again an effective approach to scattering aspects, we have exam- 
ined in a recent workJ^D] the Dirac equation with the time component of a 
vector potential in the form of a PT-symmetric square well: when the real 
depth exceeds 2m, with m the particle mass, transmission resonances at neg- 
ative energies appear as the signature of spontaneous pair production, but 
become weaker with increasing imaginary depth and negligible beyond the 
critical value at which real bound states disappear. 

In the present work , which extends the non-relativistic results of Refs.[TIj,[12j, 
we consider a scalar and vector combination of non-local separable potentials 
in the (l+l)-dimensional Dirac equation, aimed in particular at the study of 
symmetries known in their three-dimensional form as the spin and pseudo- 
spin symmetries, the latter being experimentally observed in atomic nuclei 
. Numerical results will be consistently obtained from the PT-symmetric 
version of a solvable potential originally proposed by Yamaguchi for the de- 
scription of bound and scattering states of the neutron-proton system. 

Since this kind of potential has received until now moderate attention 
within the framework of PT-symmetric quantum mechanics, and, to our 
knowledge, no attention at all in its relativistic version, we consider it worth- 
while to perform a detailed, albeit effective analysis by means of the Green 
function method described in Section 2. The scattering matrix is then stud- 
ied in Section 3 and two non-relativistic limits for the particular choices of the 
ratio of vector and scalar couplings corresponding to spin and pseudo-spin 
symmetry are discussed in Section 4. Bound states are studied in Section 5 
and numerical results obtained with a kernel corresponding to the Yamaguchi 
potential are discussed in Section 6. Section 7 is dedicated to conclusions and 
perspectives of future work. 
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2 Green function approach 

Let us start with the (l+l)-dimensional Dirac equation with a vector-plus- 
scalar non-local potential, written in units H = c — 1 



i—V(x,t) = i-ia x —+(3m)^(x,t) + (c s (3 + c v ) I dyK(x,y)^ 



+oo 



— oo 



= (a xPx + Pm)V(x,t) + (csP + c v ) / dyK(x,y)V (y,t) 

J — oo 

= (iy D *)(s,t). (2) 
A stationary wave, \& (x, t) = ^ (x) e~ lEt , satisfies the equation 

(H D V) (x) = E^> (x) . (3) 

Here, a x and (3 are 2x2 anticommuting Dirac matrices with unit square, a 2 x = 

_ ? ) = lo, which can be identified with two Pauli matrices: in the 
1/ ' 

present work we adopt the Dirac representation [2T] a x = a x = ( ^ J 

/3 = cr z = ^ q ^1 / P ar ^ cu ^ ar ^ su ited to the study of the non-relativistic 

limit of the model, cs and cy are the real strengths of the scalar potential 
and of the time component of the vector potential, respectively, with common 
"PT-symmetric kernel K(x,y) = K* (—x, —y)- 

Here, as in our previous work[20] on the one-dimensional Dirac equation 
with a PT-symmetric square well, we have the parity operator V in the Dirac 
representation 

V = e^P a z (4) 

where Pq changes x into —x and 8-p is an arbitrary phase factor. In the same 
representation, the time reversal operator T reads 

T = e ieT a z )C (5) 

where K performs complex conjugation and Or is an arbitrary phase factor. 
With the convenient choice 9? = —9-p the VT operator takes the form 

VT = P JC (6) 
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adopted also in non-relativistic quantum mechanics [11], [T2] . 

It is worthwhile to point out that formula §Z§ does not contain the most 
general Hamiltonian: for instance, we might add a pseudo-scalar interaction 
by extending the matrix of coupling strengths to cs/3 + cy + icpa x (3. The 
method of solution described in this section could be applied even to the 
most general case, but we do not consider it explicitly, because we are mainly 
interested in interaction potentials that permit decoupling the two integro- 
differential equations satisfied by the two components of so as to obtain 
a clear definition of their non-relativistic limits, as will be shown in detail in 
Section 4. 

In order to deal with a solvable model, we assume a separable kernel of 
the form 

K(x,y)=g(x)e ia *h(y)e a ", (7) 

where a and b are real numbers and the real functions g and h are even 
functions of their arguments, g(x) = g{—x) and h (y) = h(—y), so as to 
assure VT invar iance. When g = h and a = b = 0, the kernel becomes real 
symmetric and coincides with that of Ref. [22J. When g = h and a = —b the 
kernel becomes Hermitian, since in that case K{x,y) = K* (y,x). 

Now we solve Eq. ([3]) by means of the Green function method already 
used for the one-dimensional Schrodinger equation with the same type of 
potential jll] [12] - The Green function method had already been used in the 
solution of a scalar-plus- vector real non-local separable potential in Ref. [22] 
and of a pure vector potential in Ref. [23J . 

Two linearly independent Green functions, G+ (x,x') and G- (x,x'), for 
the time-independent Dirac equation are solutions to the equation 

i^~i a x-^~ + P m — {E ± ie)j G±(x, x') = S (x — x') , (8) 

where a small imaginary component e (> 0) is added to the energy, E, in order 
to remove the energy poles from the contour of the complex integral defining 
G±(x, x'), as discussed in the following part of this Section. G + (x, x') and 
G_ (x, x') are related to the Laplace transform with respect to time of the 
retarded and advanced component of the causal Green function, respectively, 
as shown in Appendix A. 

Eq. ([HD is easily solved in momentum space after introducing the Fourier 
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transforms G± (g, q') 

p+oo 

(2vr) f 

and the Fourier representation of the Dirac 5 function 



1 /* + oo /»+oo 

G ± (x, x') = —j / dqe^ / dq'e iq ' x 'G ± (g, g') (9) 

Z7T J-oo J-oo 



1 /■+00 

(J (a; -a/) = — / dqe iq{x - x ' ] . (10) 

After inserting formulae (iDT fTUl) into Eq. (jSj) we quickly obtain G-t (g, g') 
in the form 



G± (q, q') = lim e ^o+ 2vr (a x g + (3m - E =f ie) 1 5 (g + g') 
= 2tt lim e _ + ^f^t^ (g + q') 

Therefore, we obtain for the Green functions in configuration space 



G± (x, x') = G ± {x- x) = — lim / dqe iqix - xl) 

2tt ^o+_/_ oo 



1 , ; _ f ^.Ja(x-x') a xq + pm + E±ie_ 

2 ' 



q 2 + m? — {E ± ie) 

(12) 

which can be easily computed by the method of residues. Let us start with 
G + (x — x')\ after defining k 2 = E 2 — m 2 , we observe that the integrand 
has two simple poles at g x = —k — ie' and g 2 = +k + ie', where e' = eE/k. 
For x — x' > the integration contour is closed in the upper q half-plane, 
including the pole at q = g 2 , while for x — x' < the contour is closed in the 
lower g half-plane, including the pole at q = q\ with a global — sign, because 
the integration is done in the clockwise direction. The result is 

G + (x - x') = i [9 (x - x') e^*-*') (a x k + (3m + E) 

+6 {x' - x) £.-<*(*-*') (-a x k + (3m + E)] (13) 
= ±- e tk \ x - x '\ {ka x sgn (x — x') + (3m + E) , 

in agreement with Ref. [22|. In the same way we compute GL (x — x'), after 
observing that for x — x' > the integration contour in the upper q half-plane 
now includes a pole at g 3 = —k + ie' while, for x — x' < 0, the contour is 
closed in the clockwise direction in the lower g half-plane around a pole at 
g4 = +A; — ie'. The result is 

G_ (x - x') = —± [9 (x - x') e -**(*-*0 (_a s ife + /3m + £) 

+fl (x' - x) e**^ - ^ (a x Jfe + /5m + £)] (14) 
= — J- e - ik \ x ~ x I {—ka x sgn (x — x') + (3m + £7) . 



6 



Summing up 

G± (x - x') = ±-L e ± ik \*-z'\ (±ka x sgn (x - x') + (3m + E) . (15) 

It is immediate to check that 

G_ (x - x') = VTG+ (x - x') (vry 1 . 



3 Scattering matrix 

By exploiting the results of the preceding section, we can define two linearly 
independent solutions to Eq. ([3]), (x) and (x), in the implicit form 

*± {x) = *% ee (x) - L + r dx ' G ± & - + c v) IT d y K V, y) *± (y) 

= Kee ( x ) - I-Z dx ' G ± ( x ~ x ') + c v) 9 V) e iax ' dyh (y) e*»* ± (y) 

(16) 

In Eq. ( TIB]) . ^yy ee (x) is the general solution to the Dirac equation for a 
free particle, conveniently written in the matrix notation of Ref.[21j 



e ikx e -ikx \ f A±\ ( A±e ikx + B±e~ ikx 



®free( x ) \ \ ( 'A - ) ' \ B± ) \ \A±e ikx - \B±e~ ikx 

^^^^ ' (17) 

where \ = k/ (E + m) = y\E — m) / (E + m) and A± and 5± are arbitrary 
constants. It is worthwhile to point out that G± (x — x') and cs/3 + cy are 
non-commuting 2x2 matrices: therefore, their order is not arbitrary. 

After defining I± = dyh (y) e lby ^± (y), we multiply both sides of Eq. 
( Fl6l) by h (x) e lbx and integrate them over x from — oo to +oo. Remembering 
that / (q) = f_°° dxe~ iqx f (x) is the Fourier transform of / (x) and observing 
that / (x) = f (—x) implies f (q) = f (—q), we promptly obtain 

Ajx (k + b) + B±h(k - b) 
XAji {k + b)~ XBji (k - b) 
- dxh (x) e ibx j+™ dx'G ± (x - x') {c s (3 + c v ) g {x') e iax ' I± (18) 
Aji (k + b)+ B±h(k - b) 
\A±h (k + b) — \B±h (k — b) 



N± {c s (3 + c v ) I 



± 
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where N± = J_ °° dxh (x) e tbx f_ °° dx'G± (x — x') g (x 1 ) e iax ' . Therefore, the 
spinor I± is explicitly given by the relation 

once we have determined the 2x2 matrix N±, which, according to formula 
( I15p . is conveniently rewritten as 

N± = ±4 dxh (x) e ibx dx'g {x') e iax ' [ e ±ik ^~ x ') (±a x k + (3m + E) 9 (x 

+e TiHx-x') ( T a x k + (3m + E)9 (x' - x)] 

= ±± {±a x k + (3m + E) + {^a x k + (3m + E) 

(20) 

where 

/+oo f+oo 
dxh (x) e tbx / dx'g (a/) e iax ' e ±ik ^~ x '^9 (x - x') (21) 
-oo J — oo 

and 

/+oo r+oo 
dxh (x) e lbx / dx'g {x') e iax ' e^ ik{x ~ x ' ] 9 (x' - x) . 
-oo J — oo 

(22) 

It is worthwhile to point out the following symmetry relation 

N (J) _ 6> k ) = (N® (a, b, k)J (j = 1, 2) . (23) 
After introducing the linear combinations 

S± (a, b, k) = N { ± ] (a, 6, k)+N { ± ] (a, 6, jfe) , D± (a, 6, jfe) = jvj? (a, 6, A;)-ivl 2) (a, 

(24) 

with symmetry relations 

S + (-a, -6, jfe) = (S- (a, 6, jfe))* , D + (-a, -6, jfe) = (£L (a, 6, jfe))* , (25) 
iV± becomes 

iV± = ^ [Ddb^A; ± S± (/3m + E)} . (26) 
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We now specialize to the Dirac representation, already introduced in Sec- 
tion 2, a x = a x , (3 = a z . After some simple algebra, we obtain 

- ( 4 ) ' (27) 

M± = 1 2 + iV ± (c s /3 + cy) 
l±i%( Cy + Cs ) iD ± ( Cy _ Cs ) \ . (28) 
|Z>± (cy + c s ) 1 ± \XS ± (c v - c s ) J 
In order to compute I± from formula (jT9l) . we need the inverse of M± 

M -l _ _J_ ( 1 ± ( C V ~ C S ) -lD ± ( C y-C S ) \ 

M± ~ detM ± V -{D± ipv + c s ) 1 ± 1% (cv + c 5 ) J ' ^ 

with 

det M± = (1 ± | \S± (cy - c s )) (l ± |% (cy + c s )) + ^ (c 2 y - c|) 

= 1 ± 1 5± (A (<V - C S ) + i (cy + C^ + £=El (£)2 _ 5-2 ) 

(30) 

Note that, as a consequence of relations ( 12311251) . 

det M_ (-a, -6, jfe) = (det M+ (a, 6, k))* . (31) 

We are now in a position to express the asymptotic forms of the wave 
functions ty± (x) and the transmission and reflection coefficients for progres- 
sive and regressive waves in terms of known quantities. For the sake of clarity, 
let us consider ty + (x) and (x) separately. 

In order to determine the asymptotic behaviour of (x), we observe 
that 

i 

x— >±oo 



lim G+ (x - x') = L e ±ik{x - X,) ( X ^ ) ■ ( 32 ) 



Therefore, in particular 



lim x ^ +oo y + {x) = A + [ je^ + S+l _ A Je 



^j_^dx'e ik{x ~ x '^e iax 'g{x')[ J 



i 



IN. / Cy + C5 



Cy - C S 
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(33) 



If we impose the condition that \P + (x) is a progressive wave, travelling 
from left to right (L —* R), we can put A + = 1 and B + = in the preceding 
equation. After deriving from formula (fT9j) the explicit form of I + 



h(k + b) 



1 + f XS + (cy - C S ) - | \D + ( C y - C S ) 
-\D+ (cy + C S ) + A + 15+ (Cy + C 5 ) 



det M+ 

the above limit can be rewritten after some algebra in the form 



(34) 



lim^+oo ^ + (x) 



Akx 



\g (a - k)h(k + b) x 



| (c v g+c s m)+i(cy-c|)(g + -g + ) 



l+-±(c v E+c s m)- 



4 l L V C S 



(35) 



allowing us to determine the transmission coefficient, Tl^,r, since we must 
have 



lim \E' + (x) 



Akx 



(36) 



From comparison of the r.h.s. of Eqs (135^) and (I3"6"j) we obtain 

f (c v E + c s m) +i(c 2 v - c|) (5+ - £> H 



Tr^n = 1 — qf (a — k) h (k + b) — 

* l + i^(c V E + c s m) + l(4-4)(Dl-Sl) 

(37) 

In the same way we can compute the reflection coefficient, Rl^r, starting 
from 



-ikx 



lim^.oo ^ + (x) = 
dx'e^ k+a >'g (x 1 ) ( ~ 

e ikx _ ^ e -ikx~^ _|_ a 



Akx 



Cy + C S 
Cy - C S 

^ " C 5 - Cy 

Cy + C S ) A (Cy - C S ) 



(38) 



Using again formula (1341) for J + , we obtain after some simple algebra 



lim fx) 



i g {k + a) h{k + b) ( Cy + Cs 



detM 4 



A 



+ A (c S ~ Cy] 



-A 



(39) 
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where detM + is given by formula (150]) . On the other hand, we must have 
lim (x) = ( \ j e lkx + Rl-+r ( \ J e~ ifca: . (40) 



From formulae (I39p and (j40j) we promptly obtain 
R L -*r = - T g {k + a)h{k + b) 



* J l+t^(c v E + c s m) + \(c 2 v -4)(Dl-Sl) 

(41) 

In order to compare our results with those of Ref . [22] for a real symmetric 
kernel, with g (x) = h(x) = v (x) and a = b = 0, we observe that, in 
this limit, D + = and S+ = J = f^dx f*™ dx'e ik ^ x ~ x '^v (x) v (x 1 ). J 
is promptly expressed in terms of the Fourier transform of v(x), v (k) = 
J-°° dxv (x) e~ lkx . In fact 

/+oo p+oo 
dx / dx'v (x) v (x') [cos k \x — x'\ + i sin k \x — x'\] , 
-oo J — oo 

(42) 

where 

-oo , p+oo 



Jr = f dx f_ dx'v (x) v (x 1 ) cos k(x — x') 
dx f_™ dx'v (x) v (x') e lk ( x ~ x ') = v(k)v(—k) = (v(k)) 



and 



iJj = J- J R = J-{v(k)f . (44) 

Thus, in the same limit, we obtain S+ — g(k)h(k) = J — (y (k)) 2 = 
J - J R = ijj and Si - 2g(k)h(k)S + = J(J-2J R ) = -J J* = - \J\ 2 . 
Therefore, our \Tl^r\ 2 , from formula (1371) . coincides with formula (13) of 
Ref. [22J and our \Rl~^r\ 2 , from formula (jHJ, with formula (14) of the same 
reference, as expected. 

Let us now consider the second Green function, G_ (x — x'), whose asymp- 
totic behaviour is 

lim G_ (x - x ) = - l - e ^K x -') ^ ) ' ^ 
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Therefore, in particular 



lim a; _ > _ 00 ^_ (x) = A_ 



+ dx ' eik(x ~ X ' )eiaX '9 i x ') ( i \ 



i 1 



Akx 



B_ 



-A 



-ikx 



C S + Cy 

Cy - C s 



A. 



+ k (a - k) 



C S + C V A (cy - C S ) 



e lkx + 5„ 



1 

-A 
(46) 

We can impose the condition that (x) is a regressive wave, travelling 
from right to left (R —>■ L), so that 



lim m_(x)=T R ^ L 



X — > — oo 



-A 



-ikx 



(47) 



Comparison of formulae (14"6T) and (14TI) yields 



+ fflr(a-fc) „ * 



Cy - C 5 



C 5 + Cy A (Cy - C S ) 



■ L 



(48) 



In the same way 

lim^+oo ^_ (x) = A. 



+ l 2jl™ dx'e- ik{ - x - x '^e iax ' g(x') I . 



1 ' 






)_' 


" 1 




A 




— 1 


A 



e ikx + B 



1 



-ikx 



^jT/IcCC j 



A J + 1 5 (« + *0 



-A 

1 \ / cs + tv 

\ Cy - C S 

C 5 - Cy 
-(Cs + Cy) A(cy-C 5 ) 

(49) 



Since we know that 



lim (x) 

x— >+oo 



e- ikx + f?^ L 



Akx 



(50) 



we obtain 



-A 



+ y (a + fc) 



-(c S + Cy) A (Cy - C 5 ) 



1 

-A 



(51) 
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Remembering the expression of J_ from formulae (ll9ti29tl3"Ui) . and rewrit- 
ing it more compactly as 



1 - §AS_ (c y - c s ) -f ALL (tv - c s ) \ / e 



det M_ 1 ~k D - (°v + c s ) 1 - 5% (<V + c s ) 7 V AD 



, (52) 



with 6 = A_/i (Jfc + b) + (A; - b) = R R ^ L h (k + b) + T R ^ L h (k - b) and 
D =RR^Lh (k + b) — Tn^ih (k — b), we obtain from Eqs. (I48H51P a system 
of two linear equations in the unknowns & and D 

i g(a-fc). 




2detM_'^+ ~~ 

1 g(a+k) , 

2 det M- 



ty- = 1 



(53) 



with 

V± = - 1 (c 2 v - 4) is. ± d.)] e± [\ (c v - as) - J (4 - 4) ± 



(54) 
(55) 



where 

[ SS ± B ,_ |(4_4) (5_±£>_)] 

«pf } = ± [A (cy - c s ) - I (4 - 4) (£_ ± 

The transmission and reflection coefficients are thus obtained by solving 
the system ( 15*31) 



e _ s det M_ (2 det M_ + i£ (a- k)h{k+b) } +<P^ S ) 



2h(k-b) d e 

6+v ig(a-k)h(k-b) det M- (qg>-S^ 



(56) 



2h(k+b) ' d e 

with 

' »(») 1 m( e ) 



d e = 2 (det M_) 2 + ig(a-k)h(k + b) det M_ ] + <p< 

-iflf(a + fc)^(A;-6)detM_ ftp? } - <pL e) ) (57) 
(a + k)g(a-k)h(k + b)h{k- b) \ - 



Formulae ( 1561) can be further semplified by noting that 

<pf } + «pV e) = 2 Cyg + cgm - 2* (4 - 4) ivi 1} , 

jpf) _ qj(f) = _ 2 c y E+c sm + 2 - (c 2 _ c | ) ^(2) ) 
jpWfjjW) _ y^yf > = _ 2 (4 - 4) det M_ . 
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It turns out that 
"' v - -2\detM_ + i CvE t csm 



detM 



g (a - k) h (k + b) + g (a + k) h(k - b) 



+ (4 - 4) g{a-k)h{k + b) N W + g (a + k)h(k - b) N ( _l 



r(2) 



-g (a + k)g(a-k)h(k + b)h(k- b) 



2 det Mi 



(58) 

The last step is proved in detail in Appendix B. With the above result, 
formulae (151)1) are written as 



T, 



detM_ +g{a-k)h{k + b) 



;CyE+Csm 



+ K 



N. 



(i) 



R-^L 



R 



g (a — k) h(k — b) 



R-^L 



detM+ 

;^±^ + (4-4)ivi 2) 



(59) 



detM 4 

It is straightforward to verify that[12j 



T L ->r (~a, -b) 
R L ^R(~a, -b) 



T R ^ L (a,b) , 
Rr—>l { a i b) . 



The scattering matrix, S, can be defined as in Ref.pT 

S = 



Tl^r 
Rl-^r 



Rr-^l 



(60) 



(61) 



The general properties of the S matrix obtained in Ref.JTT] in case of 
V, T, or VT invariance of the Hamiltonian hold in relativistic quantum 
mechanics, too. In particular, VT symmetry of the Hamiltonian implies 



or 



s- 1 = S* 



det S\ = 1 



IT, 



L-*R\ 



IT, 



R-*L\ 



(62) 



(63) 



Im(R L ^ R R R ^ L ) = . 

Tl^r and Tr^l have the same modulus, but different phase: the latter 
property, characteristic of non-local potentials, is discussed in particular in 
Re&.[I2],[n]. 

Finally, the last of conditions fl63|) implies that Rr^l and Rl->r have the 
same phase, although they have different moduli, since unitarity is broken. 
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4 Symmetries and non-relativistic limits 



Eq. d3J) is equivalent to a pair of coupled differential equations in the two 



(m - E) (x) - ( x ) + ( Cs + c v ) dyK {x, y) ^ (y) = 

-ij^ (x) - (m + £) ^ 2 (x) + (c y - c s ) f+™ dyK (x, y) v& 2 (y) = 



For arbitrary values of the coupling strengths, cs and cy, the above equa- 
tions do not decouple; decoupling occurs when cy = ±cs- The method of 
solution described in the preceding section remains valid and the final results 
for the reflection and transmission coefficients are still given by formulae (I3T1 - 
l4Tj) for progressive waves and by formulae fl56l) for regressive waves, even if 
intermediate formulae are different. 

In 3 + 1 dimensions, the cases cy = cs and cy = — cs are examples 
of Bell-Ruegg symmetries |24j. where the Dirac Hamiltonian commutes with 
the generators of an SU (2) group, constructed with Dirac matrices and the 
momentum operator. The eigenstates of the Dirac Hamiltonian belong to 
the carrier space of the spinor representation of such a group and are thus 
doubly degenerate. When cy = cs, the members of the doublet have the 
same radial quantum number n r , the same orbital momentum I and total 
angular momentum j — I ± | (spin symmetry). When cy = — Cs, they have 
quantum numbers (n r , /, j = I + |) and (n r — 1, 1 + 2,j — I + |J, i.e. the 
same pseudo-orbital momentum 1 = 1 + 1 and pseudo-spin ? = ~, so that 
j = (pseudo-spin symmetry). The mean field of heavy nuclei exhibits an 
approximate pseudo-spin symmetry, experimentally known for many years, 
but correctly explained as a relativistic effect only few years ago[25j. At 
a phenomenological level, the approximate pseudo-spin symmetry naturally 
arises in relativistic mean field models, where the nuclear mean field is in 
practice the sum of an attractive scalar field (the a field) and of a repulsive 
vector field (the u field) of almost the same strength. At a more fundamental 
level, it can be obtained from sum rules of quantum chromodynamics in 
nuclear matter |26j. 

Let us consider the case cy = c$ = c first. We promptly obtain in this 




; in the Dirac represen- 



(64) 
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case 

(x) + 2c(m + E) f+™ dyK (x, y) ^ (y) = (E 2 - m 2 ) ^ (x) = k 2 *t 

(65) 

The above system is suited to the study of the non-relativistic limit (E — > 
m + y~ , with ^- <C m), where the first equation of system f[6"5l) . satisfied by 
becomes a Schrodinger equation with a non-local potential of strength 
s = 2c and kernel K. \I>2, being proportional to gj^i, does not obey a 
Schrodinger-like equation. 

In this limit, the transmission and reflection coefficients obtained in the 
preceding section simplify considerably. In fact, from formulae (I37H41I) we 
promptly obtain, for cy = Cs = c and E — > m + j- 



lim, k 2 T L _> R = 1 - t^ ^k-aMk+b) 

]\yt\ o R — i 2cm 9(k+a)h(k+b) ' 

11II1_ , ju2 xtr.p — — £ — ; ; . * rrn „ 

E— *m+-^ k l+ l l£01S + 

2cm 



(66) 



in agreement with formulae (153) of Ref . [1 lj . where =Tp is indicated with 
and with D + , not to be confused with the D + integral defined in 

formulae (Bill of the preceding section. It is worthwhile to recall that Ref. [TT] 
uses units 2m = 1, as is common in non-relativistic quantum mechanics. 
In the same way, we obtain, after some simple algebra 

lim T = 1 — i 2 ^ 1 g{k+a)h(k~b) 

E ^ m +h R ~" L k l+i^[-S-+g(k-a)h(k+b)+g(k+a)h(k-b)\ , . 

lim 7? = i 2cm g(k-a)h(k-b) ' v ) 

E ^ m +h k l+i^\-S„+g(k-a)h(k+b)+-g(k+a)h(k-b)\ 

which coincide with formulae (156) of Ref. [TTj , where i 2cm ^~ is indicated with 



N_, not to be confused with the N_ matrix defined in formula (J2j 

In the case cy = —cs = d ^land ^2 interchange their role, since the two 
decoupled equations now are 

-£2^2 (x) + 2d (E - m) dyK (x, y) ^ 2 (y) = (E 2 - m 2 ) ^ 2 (x) = k 2 ^ 2 

(68) 

The formulae of transmission and reflection coefficients now depend on 
E— m, to be replaced in the non-relativistic limit by the kinetic energy In 
that limit, the second equation (16"B"1) . satisfied by \l>2, becomes a Schrodinger 
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equation with an energy dependent coupling strength s (k) 
while is not solution to a Schrodinger equation. 
The final expressions are 



= c'k 2 /(2m 2 ) } 



lim 



ic'k g(k-a)h(k+b) 



lim 



2m 



c'fc g(fc+q)fc(fc+6) 
2m 1-N^S+ 



(69) 



and 




2m l+^r [-5-+g(fc-a)/i(A;+b)+s(A:+a)/i(fc-6)] 

^ c'fc g(k-a)h{k-b) ^ 

2m 1+^? [-S , _+g(fc-a)ft(fc+6)+g(fc+a)h(fc-fe)] 



(70) 



As expected, the above formulae have the same structure as those in the 
case cy = cs, with the constant strength s = 2c replaced with the energy- 
dependent strength s (k) = c'k 2 /(2m 2 ). 

For arbitrary values of cy and cs the equations AM]) do not decouple, 
unless the potential becomes local, K(x,y) = 8(x — y)V(x). As a conse- 
quence, in the particular case of a purely scalar potential, cy = 0, we do not 
obtain the pseudo-supersymmetric scheme of Ref.[13], which holds for local 
potentials only. 

Summing up, the cases cy = ±cg, reflecting the Bell-Ruegg symmetries [21] 
in one dimension, reduce the two-dimensional manifold [^1,^2] to the one- 
dimensional manifold \&i when cy = c$, or \l/ 2 when cy = — c$, the latter 
case being relevant for nuclear physics. 

5 Bound states with real energy 

The VT symmetry of the potential kernel, K, permits the general statement 
that either bound state energies are real, or that they come in complex con- 
jugate pairs: in fact, if \I/ (x) is the solution to the stationary Dirac equation 
with energy E and momentum p x 



(E - a xPx - (3m)m (x) - (c s (3 + cy) / dyK {x, y)V (y) = , (71) 




—00 
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VT^fr (x) = v &* (—x) is solution to the Dirac equation with energy E* and 
momentum — p* x 

+00 

(E* + a x p* x - {3m)VT^ (x) - (c s f3 + c v ) j dyK*(-x, -y)VT^ (y) = . 

—00 

(72) 

We now treat in particular bound states ^bs (x) with real energy and 
imaginary momentum p x = —p* x and investigate the relation between ^bs (x) 
and VT^bs (%)■ From now on, the quantum number k is no more real and 
positive, as defined in Section 2, but complex. 

The Green function formalism permits not only derivation of scattering, 
but also of bound state wave functions. As is known, bound state energies 
are located in the interval — m < E < +m, where the square of the mo- 
mentum, k 2 = E 2 — m 2 , is negative, i. e., k = ik is imaginary. Bound 
state wave functions can be obtained by analytic continuation of one of the 
two independent scattering solutions, e. g. (x) from formula f|T6l) . to the 
positive imaginary k axis, 1. e. we can take k = \Jm? — E 2 > 0, and impose 
the boundary conditions lim^-toc, \&+ (x) = 0. Owing to the fact that the 
Green function G + vanishes at x = ±00 when k = ik, we must get rid of the 
free-wave contribution, by putting A + = B + = 0. We thus obtain 

+00 

* 6s (x) = - I dx'g (x') e iax 'G + (x - x') (c s a z + c v ) I + . (k = ik) (73) 



Remembering expression (|T5l) for G + (x — x') , (%) can be put in the 
form 

®bs (?) = -± (e-*»Ji (x) (ikcr x + ma z + E) 

_ 1 1 (74) 

+e kx I 2 (x) {-ika x + ma z + E) j (c s o z + c v ) I + , 

where 

Ji(x)= r dx'g(x')e( ia+I >' (75) 



and 

/•+00 _ 
J 2 (x) = / dx'g (x f ) e ( ia - k y . (76) 



The normalizability of ^bs is easily cheched by noting that lim^-toc, e kx X\ [x] 
lim^ioo e kx l 2 (x) = 0. 
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From definitions (l75H76p . remembering that g (x r ) = g {—x'), it is easy to 
verify that 

J 2 (x) = XI (-x) = m (x) . (77) 

The integral equation f!73|) allows us to compute bound state energies, 
too. By multiplying both sides by h (x) e tbx and integrating them over x 
from — oo to +00, we obtain 

/+00 r+00 
dxh (x) e ibx / dx'g {x') e iax 'G + {x - x') (c s a z + c v ) 1+ , 
-00 J —00 

(78) 

or, remembering definition (1281) of matrix M + 





+00 

dxh(x)e 




■00 


( M + 




[ M 21 


M 22 J I 



ibx 



/+00 
dx'g{x')e mx 'G + {x - x')(c s a z + c v ))I+ 
■00 

= . (79) 



Note that, since M± (ik) = VTM± (ik) (VT) \ if I± is solution of Eq. 
9]), VTI± is solution, too. 
The necessary condition for a non-trivial solution of the above equation 

detM + = 1 + c v E + c s m S+ + Ml^ll (D% - Si) = , (80) 
Vm 2 - E 2 4 V + +; v ) 

where S+ and D + are functions of k(E), fixes bound state energies as the 
roots of the equation in the interval — m < E < +m. Not surprisingly, bound 
states correspond to poles of the transmission coefficient T L ^ R ( 157|) . 

Eq. f!79l) allows one to express the ratio of the components of I + in 
terms of M + matrix elements. In general, one observes that M^ 2 = 1 + 
i| (cy — cs) S + 7^ 0, so that one can exploit the second eq. fl79l) . which gives 
72 = _ (Mf/Mf) I\ and ^ bs (x) can be written as 

il f _-kx-r 1 \ ( E + m ik 



+e ^ Mx) ( E +™ F ~ l ~ k )){ CV T S c\ ){ ■ ^ 

(81) 

In the above formula, the modulus of 1} is to be determined from nor- 
malization of the wave function ty^s It is easy to check that ty^s (#) is 
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eigenstate of VT, since the matrices in curly brackets are PT-symmetric, 
owing to relation (1771) and the ratio r = — ~ - ^Y^^^t^ is real. In fact, 

from definitions (12111221), we see that, for k = ik, ivj 2) (ik) = (n^ (ik)\ ; 

thus, S + (ik) = (ik) + ivj 2) (ik) = 2Re (n^ (ifc)) is real, D + (ik) = 

N { + ] (ik) - ivj 2) (ik) = 2ilm (n£ ] (ik)^j is imaginary and A = ik/ (E + m) 

is imaginary, too, so that, as a final result, r is real and VT^bs (x) = ^bs (x), 
if l\ is chosen to be real. 

Considering that D + and S + do not depend on cy, or C5, Eq. (IHUl) can 
also be used to determine either potential strength (cy or c^), provided the 
other is fixed, in particular set to zero, in such a way to obtain a bound state 
at a given energy E in the (— m, +m) range. In this procedure, however, 
PT symmetry is not automatically preserved, since Eq. f[5U|) is of second 
degree in the unknown potential strength and might have a pair of complex 
conjugate solutions. 

We have derived our expressions for bound-state wave functions starting 
from ty + (x), but we could, alternatively, start from (x) and determine 
the constants A_ and B_ from the boundary conditions linLj^-i-oo (x) = 0. 
In this case, both A_ and B_ must be different from zero because of the 
asymptotic behaviour of G_ and appear as the solution of a system of two 
homogeneous linear equations. The condition for a non-trivial solution of the 
system yields again the equation detM + (ik) = 0, as expected, with detM + 
written in terms of det M_ according to formula (lo'Sj) . 



6 The Yamaguchi potential 

As an example of application of the formalism developed in the preceding 
sections, we now work out in detail a one-dimensional PT-symmetric version 
of the Yamaguchi potential[27|, originally aimed at describing bound and 
scattering states of the neutron-proton system. We assume 

g (x) = exp (— c \x\) , h (y) = exp (— d \y\) , (—00 < x,y < +00) (82) 

with c and d positive constants, so that the Fourier transforms are 

( g ) = A 2 1„2 ' ^ = rl2 1 „n ■ (~°° < g» < +°°) ( 83 ) 
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and the PT-symmetric kernel reads 



K(x,y) = e -4x\+iax e -d\y\-Hby^ 



(84) 



with a and b real constants. With the above definitions the basic integrals 
N±^ from formula (12T1) and N±^ from formula ( 1221) . as well as their linear 
combinations S± = N± + N± and -D-t = Af±^ — AT^ , can be computed by 
elementary methods. We only quote the final results 



N^ = --g(a T k) 



(a + b)c+{c + d){aTk) g(a^k)h{b±k) ( .b±k 

d 

(85) 
.b =F k 



(a + bf + (c + d) 2 



AT {2) _ i~ ( u ,(a + b)c+(c + d)(a± k) g (a ±k)h{b^k) 

' ~~ ~9 ^ ' ' » \2 , / _ — J\2 ' 2 



(a + &r + (c + d)' 



d 

(86) 

Formulae (186]) become particularly simple when applied to the analysis 
of bound states: in this case, we already know from the previous section 

that ivf > (ik) 
and D+ = N_ 



N±^ (ik) ) . Therefore, Sj, 



2Re N 



r(l) 



(1) 



AT 



(2) 



2ilm {N^h the left-hand-side of Eq. (ED]) thus 
becomes real in the interval — m < E < +m: bound state energies are roots 
of the real equation 



detM 4 



1 + 



2 (c v E + c s m) 
\/m 2 — E 2 



Re - (4 - 4) 



N 



(i) 







(87) 



If we put cs = in the above equation, this allows us to derive the 
strength cy at which the purely vector potential has a bound state at given 
real energy E; in fact, Eq. (157)1 can be considered as a quadratic equation in 
cy, with real solutions 



ERel N 



± 



m 2 —E 2 



+ 



iV 



(1) 



AT 



(i) 



(88) 



Note that one of the two solutions for cy is always positive. 
In order to complete the discussion of bound state wave functions, we 
give the corresponding expressions of integrals (1751) and [761) . obtained by 
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elementary integration: 

e-^ (*) = 9 (-*) gg£ + 61 



c—k+ia 



c+k+ia 



c+k—ia 



kx ^ e ( — c+ia)x g~kx 



-c+k+ia 



+ e(x) 



e (—c+ia)x 
c+k—ia 



(89) 



While possible bound state wave functions with real energy are eigenstates 
of VT , scattering wave functions never are, but show some interesting pe- 
culiarities related to transmission resonances when cy = ±cs, which makes 
it worthwhile to focus our numerical analysis on that cases. Figure 1 shows 
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Figure 1: Square moduli of transmission and reflection coefficients vs. energy 
(in units of particle mass) for Yamaguchi potentials with a = 2, b = 1, 
c — d — 1 and cy = Cs = 5m (upper panels), or cy = — Cs = 5m (lower 
panels) . 



the square moduli of transmission coefficients, \Tl^r\ 



: |T^ L | 2 EE |T| 2 , 

and of reflection coefficients, \Rl~*r\ and \Rr^l\ , as functions of total en- 
ergy E for the following choices of potential parameters: a = 2, b = 1, 
c — d — 1 and cy = Cs = 5m (upper panels), or cy = —Cs = 5m (lower 
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panels). E ranges from —5m to +5m, but the coefficients are not calcu- 
lated in the — m < E < +m interval, where they might have poles corre- 
sponding to bound states. When cy = cs, there is a sharp transmission 
resonance at E — —in, which appears at E — +m when cy = —cs, as ex- 
pected from the relation Tl->r (cy, cs, E, k) = Tl^r (— cy, cs, —E, k). These 
zero-energy resonances, called half-bound states, have T L ^ R = T R _> L = 1 
and Rl^r = Rr^l = 0. In both cases, the reflection coefficients show the 
handedness discussed in Ref.[T0]: the potentials behave as absorptive for pro- 
gressive waves (\Tl-+r\ 2 + |-Rl^_r| 2 < 1) and generative for regressive waves 
( \Tr^l\ 2 + \Rr-^l\ 2 > 1 )• This pattern depends on the (common) sign of a 
and b : in fact, owing to the form (J7J) of the kernel, where g(x) = e~ c ' x ' and 
h(y) = e~ d \ y \ are even functions of their arguments, changing a into —a and 
b into —b is equivalent to a parity transformation ( x — > —x and y — > — y ), 
namely 

T l ^r (~a, -b) = Tr^, l (a, b) , , gQ , 
Rl^r (~a, -b) = Rr^l (a, b) . 

In our case, with a = —2 and b = — 1, the potential would become gener- 
ative for progressive waves and absorptive for regressive ones. Handedness, 
however, is not a general rule: Figure 2 shows transmission and reflection 
coefficients for a = —2, b = +1, c = d = 1 and cy = cs = 2m (upper panels), 
or cy = —cs = 2m (lower panels). In this case, both \T L ^ R \ 2 + |_R^/j| 2 and 

1 2 i 1 2 

ITij_i.il + \Rr-*l\ may be ^ 1 in different energy intervals. 

As for bound states, they exist only in the lower panel cases of Figs. 1,2: 
when cy = —cs = 5m, a = 2, b=l,c = d= l( Fig.l ) there is a real bound 
state with energy = +0. 3835m, when cy = —cs = 2m, a = —2, 6=1, 
c = d = 1 ( Fig.2 ), there is a real bound state at = 0.1815m. If a and b 
change, the bound states change their energies, but they do not disappear, 
unless |a|, |6| — > +oo. In this latter case, the kernel K(x,y) undergoes such 
rapid oscillations in the x, or y directions that it becomes negligible on the 
average and cannot sustain bound states any more. In this limit, |T| — > 1 
and \R\ -> 0. 

As far as bound states are concerned, the structure of Eq. (157)1 shows 
that only when cy = det M + does not depend on E, but on k only, so 

that, if kbs is a solution of detM + (k) = 0, both energies = ± \J m? — k\ s 

are acceptable. This is shown in Fig.3, where det M + Ck) = is solved 
graphically for a scalar well of strength cs = —m, c = d = 1 and various 
values of a = b. With increasing the latter phases, the two bound state 
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Figure 2: Square moduli of transmission and reflection coefficients vs. energy 
(in units of particle mass) for Yamaguchi potentials with a = —2, 6=1, 
c = d = 1 and cy = cs = 2m (upper panels), or cy = —cs = 2m (lower 
panels). 



energies quickly tend to the thresholds of continuum, ±m. For instance, 
when a = b = 10, e& s = ±0. 999999923m and, in the continuum of scattering 
states, the potential is almost reflectionless. 



7 Conclusions and perspectives 

In this work we have studied non-local PT-symmetric potentials in the one- 
dimensional Dirac equation. Owing to the fact that the definition of the S 
matrix adopted in our previous work [11] dedicated to non-relativistic quan- 
tum mechanics is valid also in the relativistic case (see e.g. Ref.[28]), we have 
used in the present work general properties of the S matrix under V, T and 
VT transformations derived in Ref.[TT]. There are, of course, kinematical 
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Figure 3: Denominator of the transmission coefficient for a scalar well with 
Cs = —m, c = d = 1 at various values of a = b. 

differences between Schrodinger and Dirac formulations: in the latter case, 
total energies E can be either positive or negative; scattering states have ei- 
ther E/m < —1 or E/m > +1, while bound states are found in the interval 
-f < E/m < +1. 

The separable potential we have studied is very flexible, since, for in- 
stance, it permits determining the real vector strength cy ( with scalar 
strength c$ = 0) that yields a bound state at an energy E arbitrarily chosen 
in the [— m, +m] interval (see Eqs. (l8UH8Tj) ). 

Moreover, starting from the real kernel with real coupling strengths Cy 
and cs and a = b = 0, one can extend it in a natural way to the generalized 
Hermitian case, with g = h and a = —b, and, finally, to the PT-symmetric 
case, with g and h even functions of their arguments and arbitrary a and b. 

The specific choice of form factors g(x) = e~ c ' a; ' and h(y) = e~ d ' y ' yields 
in the non-relativistic case transmission and reflection coefficients that are 
rational functions of momentum, k, since they can be written as ratios of 
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polynomials in k. This opens the way to an algebraic search for zeros of 
denominators, providing information on bound states, and of numerators, 
e. g. of reflection coefficients ( transparency at given momentum k), or 
transmission coefficients ( total reflectiveness at given k). 

In the relativistic case the functional dependence is more involved, due 
to the square root dependence on k of energy E = \Jk 2 + m 2 . Nevertheless, 
it is interesting to remark that, in addition to the study of properties of T 
and R at given cy and cs, one can study specific properties like absence of 
reflection or of transmission at given k as functions of cy and cs'- this can 
be easily done since transmission and reflection coefficients are, respectively, 
second order polynomial in cy and/or cs over second order polynomial and 
first order over second order. 

Study of the zeros of the denominators has already been mentioned in 
connection with bound states. In the present work, we have made an effective 
approach to VT symmetry, allowing for unitarity breaking of the scattering 
matrix. The search for a Hermitian equivalent description would imply the 
definition of a charge conjugation operator C, in the spirit of Ref. [3], or a 
metric operator 77+, according to Ref. [6] and the study would be far from 
trivial. To our knowledge, 77+ in relativistic problems involving scattering 
states has been exactly determined until now only for a non-Hermitian form 
of the Klein-Gordon equation, either free[29j, or with a minimally coupled 
electromagnetic field [5U]. 

This kind of more fundamental study, however, would be more appropri- 
ate to finite-range potentials with exact VT symmetry, i.e. having a purely 
real discrete spectrum with eigen-f unctions that are eigenstates of VT and 
reflect ionless in the continuum, as discussed in ref.|llj. This could not be 
pursued for non-local potentials, but it could work for the "PT-symmetric 
generalization of local scalar, or pseudoscalar reflect ionless potentials, like 
those constructed in Refs.[3l]-[32]. 



Here we explain the connection between the time-independent Green func- 
tions used in the present work and the time dependent ones, which are solu- 
tions to the equation 



—ioi x — — h (3m — i— I G (x, t; x', t') — 5 (x — x') 5 (t — t') . (91) 
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We know from textbooks [22] that particular solutions to Eq. (1911) are the 
retarded component of the causal Green function 

G c ret {x, t; x, t') = 9(t- t') G+ (x - x , t - t') (92) 

and the advanced component 

G c a dv. (x, t; x', f) = 9(t'- t) G_ (x -x',t- , (93) 

where 9 (r) = 1 for r > and 9 (r) = for r < 0. 

By inserting formulae (I92H93P into Eq. (|9T|) . we obtain 



9 (± (t - f )) (-i^l + (3m - zf ) G± (x - x', t - t') 
+i5 (t - *') G± (x -x',t- t') =8(x-x')5(t- t') , 



(94) 



where we have exploited the well-known relation (± (t — t')) = 5 (t — t'). 
Let us multiply both sides of Eq. f[9~4|) by exp (iS (t — t')) and integrate 
them over u = t — t' from -oo to +oo, with £ a complex number whose 
imaginary part is chosen in such a way that the integral exists: we must 
assume £ + =E + it for G + and £-=E — it for with e > 0. Thus, G + 
satisfies the equation 

( d d \ 

due l8+u I -ia x — + (3m - i— J G + (x — x', u)-iG + (x - x', 0) = 5 (x - x') . 
\ ox ou J 

(95) 

The above equation can be simplified by integrating by parts the third 
integral on the left-hand side 

r+oo o 

-i / due i£+u — G+(x-x',u) = -i|e <£ + u G+(a;-ar',u)|^' 00 

r+oo r+oo 

-£+ due ie+u G + (x-x',u) = %G + {x - x', 0) - £ + / due i£+u G + {x - x', u) 
Jo Jo 

and the function 

r+oo r+oo 

G+ (x -x')= / due lS+u G + (x - x', u) = due iEu ~ eu G + (x - x' , u) , 
Jo Jo 

(96) 

which is nothing but the Laplace transform of G + (x — x',t — t') with respect 
to time, satisfies the equation 

\ iax ~Hx + I 3m -( E + ie )) G + ( x -x') = 5(x- x') (97) 
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and can be identified with the Green function corresponding to the complex 
energy E + ie for the time-independent Dirac equation. 

We can proceed in the same way for G_ (x — x',t — t'), after introducing 
the complex energy £_ = E — ie 



/ due t£ ~ u ( — ia x — — h dm — i-— I G- (x — x' u)+iG- (x — x', 0) = 5 (x — x) . 
J-oo \ 9x duj 

(98) 

After integrating by parts the third integral on the l.h.s. of the above 
equation and defining the Laplace transform with respect to time of G_ (x — x',t' — 

G_ (a; - x) = / due i£ - u G. (x - x', u) = di;e* a G- (x - x', -v) , 



(99) 

we arrive at the equation satisfied by the time-independent Green function 
GL(x-x') 

d \ 
-ia x — +/3m-(E- ie) GL (x - x) = 5 (x - x) . (100) 
ox J 

Eqs. (IDTI fTUOj) coincide with Eq. (JE} of the text. 



B Appendix 

Formula (158]) can be easily proved from definition (1301) . according to which 
det M+ - det M_ = i ^ E ^m + S ^ + 4_4 (#2 _ D 2 _ 32 + 32 j 

= l£x E±csrn ^(D + N (2) + ^(1) + ^ + ^ _ ^ ^(1)^(2) _ rfX jfMj > 

(101) 

where integrals N± and N± are defined by formulae (1211) and ( 1221) re- 
spectively. We promptly obtain from the definitions 

iVj 1} + iVj 2) + N ( y + iVi 2) = 2 dxh (x) e ibx dx' 5 (x') e iaa; ' cos (A; (x - a/)) 9 
+2 J"_ 00 <ix/i (x) e tbx dx'g (x') e mx ' cos {k (x — x')) 9 (x 7 — x) 
= dxh (x) e ibx dx'g (x') e ias ' ( e ifc (^') + e -*(*-*')) 
= g(a-k)h(b + k)+g(a + k)h(b-k) . 

(102) 
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and 

N W N m _ jvj 1 ) aH 2 ) = -g ( a _ k) h (b + k) g (a + k) h (b - k) , . 
+g{a-k)h{b + k)N W +g{a + k)h{b-k)N m . [ ' 

In deriving the last expression, the relation 9 (— x) — 1 — 6 {x) has been 
used in the integrands. Inserting the right-hand sides of formulae (1102111031) 
into formula (11 01 1) yields formula (158]) of the text. 

It is worthwhile to point out that the definitions we have used and, conse- 
quently, the relation between det M + and det M_ are also valid for complex 
k, in particular for k = ik, with k > 0, characterizing bound states with real 
energy. 
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